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I. INTRODUCTION 

Composed primarily of hydrogen (91% by number) and 
helium (9%), with trace amounts (0.1%) of heavier ele- 
ments, the interstellar medium plays a vital role in the 
cycle of stellar birth and death and galactic evolution. 
Not only do the properties of the interstellar medium 
govern the formation of new stars, but through their ra- 
diation and the matter and kinetic energy from their out- 
flows and supernovae, these stars in turn determine the 
properties of the interstellar medium from which the next 
generation of stars will be born. One of these feedback 
processes, the subject of this review, is the large-scale 
ionization of the medium by the youngest and most lu- 
minous stars, the O stars. Even though they are located 
near the galactic midplane in rare, isolated regions of 
star formation and often surrounded by opaque clouds of 
neutral hydrogen, the Lyman continuum radiation from 
these hot stars is somehow able to propagate large dis- 
tances through the disk and into the galaxy's halo to 
produce extensive ionization of the interstellar hydrogen. 
The study of this wide-spread plasma has impacted our 
understanding of the dynamic interstellar processes oc- 
curring in galaxies. 

This area of study began more than four decades ago, 
when |Hoyle and Ellis (1963) proposed to a skeptical as- 
tronomical community the existence of an extensive layer 



of warm (10 K), low density (10~ cm~ ) ionized hy- 
drogen surrounding the plane of our Galaxy and having 
a power requirement comparable to the ionizing luminos- 
ity of the Galaxy's O and B stars. Their conclusion was 
based upon their discovery of a free- free absorption signa- 
ture in the observations of the Galactic synchrotron back- 
ground at frequencies between 1.5 and 10 MHz carried 
out by radio astronomy pione ers Grote Reber and G. R. 
A. Ellis at Hobart, Tasmania ( ^Ellis et al. , 1962 Reber and 



Ellis 



1956|. The idea that a significant fraction of the 



Galaxy's ionizing photons, produced primarily by rare, 
massive stars residing near the Galactic midplane, trav- 
eled hundreds of parsecs throughout the disk to produce 
wide-spread ionization conflicted with the traditional pic- 
ture in which the neutral atomic hydrogen (the primary 
component of the medium by mass and opaque to hydro- 
gen ionizing radiation) fllled much of the interstellar vol- 
ume. However, less than a decade later, the dispersion 



of radio signals from newly discovered pulsars (Guelin 



1974 Hewish et al. 1968 ) plus the detection of faint op- 



warm ionized hydrogen as a major, wide-spread compo- 
nent of our Galaxy's interstellar medium. Two decades 
later, deep Ha imaging with CCD detectors revealed that 
similar warm plasmas also permeate the disks and halos 



of other galaxies (Dettmar 1990; Rand et al. 1990). 



The large mass, thickness, and power associated with 
this ionized layer has modifled our understanding of the 
composition and structure of the interstellar medium and 
the distribution of Lyman continuum radiation within 
galaxies. Its weight is a major source of in t erstel lar pres- 
sure at the midplane ( [Boulares and Cox , 1990), and it 
could be the dominant state of the inters t ellar m edium 
1000 pc above the midplane (Reynolds 1991a). An 



accurate understanding of interstellar matter and pro- 
cesses thus requires as thorough a knowledge of this dif- 
fuse ionized gas as the other principal components of 
the medium. For a general survey of the interstellar 
medium, we refer the reader to The Interstellar Environ- 
ment of Our Galaxy by Katia Ferriere (2001) and The 
Three Pha se Interstellar Medium Revisited by Donald 
Cox ( 2005 ) . Brief review articles on a variety of current 



interstellar medium topics can be found in How Does the 
Galaxy Wo rk?, a conference proceedings edited by jAlfaro 
et a/. I (12004). 



Although the nature of this low density plasma is not 
yet fully understood, signiflcant progress has been made 
in characterizing its properties and the source of its ion- 
ization. In ^nj we present a brief account of observa- 
tions of the warm ionized medium (WIM) in our Galaxy, 
the Milky Way, including results on the physical condi- 
tions in the WIM. In ^III and ^IV, we review progress on 



understanding this gas and its close connection to star 
formation activity and O stars through studies of other 
galaxies. O star photoionization models in a clumpy in- 
terstellar medium are discussed in with a discussion 
of supplemental radiation from hot-cool gas interfaces in 



^ VI A list of some unanswered questions and future chal- 



tical emission lines from the diffuse interstellar medium 



(Reynolds, 1971 Reynolds et al. 1973) flrmly established 



lenges is presented in ^ VII Throughout, we keep the con- 
vention that WIM refers to the warm ionized medium in 
our Galaxy and DIG refers to this diffuse ionized gas 
in other galaxies. Also, we use H+ to refer to the ion- 
ized hydrogen in the low density, wide-spread WIM/DIG 
and H ll to denote the gas in the localized, higher den- 
sity, "classical H ll regions" immediately surrounding hot 
stars. For reference, a list of acronyms and terms with 
explanations is provided at the end of this paper. 



II. VIEW FROM THE INSIDE: DIVING INTO THE 
PHYSICS OF THE WARM IONIZED MEDIUM IN OUR 
GALAXY 

A. Basic characteristics of the WIM and diagnostic tools 

Our location within the disk of the Galaxy provides an 
opportunity to explore close up and in detail the distri- 
bution and physical properties of this ionized medium, 
including its ionization state and temperature. The ba- 
sic features of the WIM are not very different from those 
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first proposed by Hoyle and Ellis (1963). Temperatures 



range from about 6000 K to 10 000 K, and in the solar 
neighborhood, its average hydrogen ionization rate is ap- 
proximately 4 X 10^ s~^ within a one cm^ column extend- 
ing perpendicular through the Galactic disk, abo ut l/8th 
that available from ste l lar ion izing photons (e.g., |Madsen| 
et al. , 2006 Reynolds, 1984). The amount of ionization 
increases toward the Galactic center (e.g., Madse n and] 
Reynolds! |2005). In other galaxies, the average ioniza- 



tion rate is observed to be about 1/2 that available from 
the stars (see ^IV). 



Two fundamental parameters of ionized gas along 
any line of sight, 5, are the dispersion measures (DM 
= / Ueds)^ derived from pulsar observations, and the 
emission measure (EM = J nenii+ ds ^ J nlds), de- 
rived from the intensity of the hydrogen Balmer-alpha 
(Ha) recombination line or from the amount of free-free 
(Bremsstrahlung) emission or absorption. Comparison of 
these measurements along common lines of sight indicate 
that the H+ is clumped into regions having an average 
electron density, Ug = 0.03-0.08 cm~^, and filling a frac- 
tion, / ^ 0.4-0.2, of the volume within a 2000-3000 pc 



thick layer about the Galactic midplane ( Hill et al 2008 
Reynolds, 1991b). Data also suggest that the filling frac 
tion increases from / ^ 0.1 at the midplane to / > 0.3- 
0.4 at \z\ = 1000 pc (Berkhuijsen et al 2006| iGaensler 



et al] 2008 Kulkarni and Heiles, 1987 Reynolds, 1991b). 



The large, 1000-1800 pc scale height, significantly larger 
than that of the neutral hydrogen layer, has been de- 



duced from both pulsar observations (Gaensler et al 



2008 Reynolds 1989) and from the rate of decrease in 



the Ha intensity with increasing Galactic latitude for the 



gas associate d with the Perseus spiral arm (e.g., Haffner 
1999| ). The WIM accounts for 90% or more of 



et al 



the ionized hydrogen within the interstellar medium, and 
along lines of sight at high Galactic latitude (i.e., away 
from the Galactic midplane), the column density of the 
H+ is approximat ely 1/3 that of the neutral hydrogen 
( |Reynolds|[l991aD . 

Although originally detected by radio observations, 
subsequent developments in high-throughput Fabry- 
Perot spectroscopy and CCD imaging techniques have 
demonstrated that the primary source of information 
about the distribution, kinematics, and other physical 
properties of the WIM is through the detection and study 
of faint interstellar emission lines at optical wavelengths. 
For example, the distribution of the H+ is revealed by 
its interstellar Ha (A6563) recombination line emission, 
which covers the sky. Several deep Ha surveys have given 
us our first detailed view of the distribution and kinemat- 



i cs of this gas. Through CCD imaging, |Dennison et al 
(11998') has provided partial coverage of the total Ha in- 



tensity in the northern sky, while Gaustad et al. (2001) 



covered the southern sky at arcminute resolution and sen- 
sitivities around 1 R (10^/47r photons cm~^ s~^ sr~^; 
EM = 2.25 cm-^ pc at 8000 K). The Wisconsin H- Alpha 
Mapper (WHAM) carried out a velocity resolved survey 
of the northern sky (dec. > 30°) at 12 km s~^ spectral 




FIG. 1 Composite all-sky velocity-integrated Ha map. Data 
from VTSS, SHASSA, and the WHAM northern sky survey 
(WHAM-NSS) have been combined to produce these very 
deep (EM > 1 cm~^ pc) emission maps. The two Aitoff- 
Hammer projections are centered at (top) i = 0°and (bot- 
tom) £ = 180°. Areas covered by the imaging surveys have 
arc-minute resolution while those only surveyed by WHAM 
have one-degree resolution. Emission is predominantly from 
the Galaxy, although the imaging surveys can contain other 
bright, extended sources from the local universe (I i'lsrI S, 500 



km s-^ LMC, SMC, M31, etc. 
I2QQ31 



Adopted from 



Finkbeiner 



resolution, 1° spatial resolution, and ^ 0.1 R sensitivity 
(iHaffner et al] [20031 [Tuftel \m7\. The WHAM maps 



show WIM emission in virtually every beam, with faint 
loops, filaments, and blobs of emission superposed on a 
more diffuse background. Finkbeiner| ( |2003[ ) has com- 



bined these three surveys to form a composite all-sky 
view of the velocity-integrated Ha (Fig. [T]). 

With these new maps and methods for detecting faint 
emission lines, we are now able to investigate the physical 
conditions of the WIM, its relationship to other compo- 
nents of the interstellar medium, and to sources of ion- 
ization and heating within the Galaxy. In particular, 
standard nebular line diagnostic techniques can now be 
employed to examine the physical conditions in the gas. 
In the low density (~ 10~^ cm~^) environment of the 
WIM, the collisional excitation of an ion to a metastable 
state 2-3 eV above ground by the thermal (~ 10^ K) 
electrons is followed by the decay back to the ground 
state via a "forbidden" optical transition. Specifically, 
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the ion's excitation rate oc rie T~^-^ ex.p{—E/kTe)^ 
where rii and rie are the volume densities of the ions 
and electrons, respectively, Tg is the electron tempera- 
ture, and E is the energy of the metastable state above 
ground. Because thermal equilibrium between electrons 
and ions is very rapid, the temperature of the ions Ti = 
Te (e.g., |Spitzerl|1978 V Thus a variation in the photon 
emissivity of a forbidden line from one direction to the 
next traces variations in the temperature, density, and 
abundance of the ion. The effects of density variations 
can be eliminated by dividing the forbidden line inten- 
sity by the H-recombination line intensity, both of which 
are proportional to the product ng. From the intensi- 
ties of lines from a number of different ions and atoms, 
it has been possible to study separately variations in the 
temperature and the ionization state within the emitting 
gas. For many years, these diagnostic techniques have 
been applied to a variety of astrophysical plasmas (see. 



e.g.,|Davidson and Netzerl|1979 


|Dopita and Sutherland , 


20031 Ferland 


, |2003t lOsterbrod 


kl |1989t lOsterbrock and 


Ferland,|2006 


for in depth discussions), but only more 



recently has it been possible to use them to study the 
much fainter WIM emission lines. 

For example, in the WIM, the forbidden lines [S ll] 
A6716 and [N ll] A6584 are found to have intensities with 
respect to Ha that range from a few tenths to unity or 
higher, significantly larger than what is observed for the 
bright, classical emission nebulae (i.e., H ll regions) im- 
mediately surrounding O stars. This implies that the 
physical conditions in the WIM differ significantly from 
conditions in classical H ll regions. In addition, because 
their intensities are comparable to Ha, it has been pos- 
sible to map these lines over large parts of the sky (e.g.. 



Madsen etaL\\2006\ . Other lines, such as [N ll] A5755, 
He I A5876, [Oln] A5007, and [O i] A6300, are much 
fainter and have been studied only in a few select di- 
rections. These observations have helped to characterize 
the ionization and temperature of the WIM as well as 
other extended ionized regions of the Galaxy. Results 
reveal that not only are the temperature and ionization 
conditions of the WIM significantly different from the 
conditions in classical O star H ll regions, but that the 
conditions within the WIM itself vary considerably from 
one direction to the next and even along a single line of 



sight (Madsen et a/., 2006). 



B. Ionization state 

The strength of the ionizing radiation field responsible 
for the WIM can be probed by measuring the Ha sur- 
face brightness of neutral hydrogen (H l) clouds and by 
measurin g the hydroge n ionization fraction H+/H within 
the WIM. [Field] (|1975|) pointed out that an H I cloud im- 



mersed in an ionizing radiation field will have a skin of H+ 
with an emission measure that is directly proportional to 



47r J 2 X 10^ photons cm~^ s~^ could account for most 
of the WIM's ionization. This flux implies an ionizing 
photon density to electron density ratio (the ionization 
parameter) of 10~^ to 10~^, which is one to two orders 
of magnitude smaller than the ionization parameter in 
classical O star H ll regions. However, values in this 
range still imply that the hydrogen is nearly fully ionized 
within the WIM. This is confirmed by more direct mea- 
surements of the hydrogen ionization fraction from the 
detection of neutral oxygen emission. 

In theory, directly measuring the degree of H-ionization 
within warm, ionized gas is simply a matter of observing 
the [O l] A6300 emission line, which is produced by colli- 
sions of neutral oxygen with thermal electrons within the 
WIM. The first ionization potential of O is quite close to 
that of H (13.595 eV and 13.614 eV, respectively) and 
the large H+ + ^ H^ + 0+ charge-exchange cross 
section keeps 0+/0 nearly equal to H+/H. Electron 
energies in Tg ~ 10^ K gas are sufficient to excite the 
~ 2 eV (^P — ^D) transition that results in the [O l] 
A6300 emission. Therefore, the intensity of this line rel- 
ative to Ha is directly related to the amount of O^, and 
thus H^, relative to H+ in warm ionized gas (Reynolds 



et al, 1998). In practice, nature conspires to make this 
observation very difficult, because [O l] is also one of the 
brightest terrestrial emission lines in the night sky. Nev- 
ertheless, high-sensitivity, high-resolution spectroscopic 
measurements with WHAM have managed to resolve the 
Galactic emission from atmospheric emission to provide 
reliable measurements in a few select directions. These 
observations indicate that H+/H > 90% fo r T > 8000 K 
(iHausen et a/.l 120021 [Reynolds et a/.l 119981). 



The time scale for recombination at a typical WIM 
density of 0.1 cm~^ is ^i:^ 1 Myr. This is shorter than the 
lifetimes of O stars, the presumed ionizing sources, which 
implies that the photoionization rate of the neutral hy- 
drogen atoms within the WIM is roughly balanced by the 
rate of hydrogen recombination. In this case, the limit on 
H+/H implies an ionizing flux > 10^ photons cm~^ s~^, 
consistent with the ionizing photon flux derived from the 
Ha surface brightness of H I clouds. 

Regarding heavier ions, observations reveal that in the 
WIM ions are generally in lower states of ionization than 
in classical O star H ll regions (e.g., Madsen et al. , 2006). 
The reason is not yet fully understood, since Lyman con- 
tinuum photons emitted by massive O stars are almost 
certainly the primary source of ionization for the WIM 
(see ^III and ^IV| below). The lower ionization state could 
be due to a softening of portions of the spectrum as the 
radiation travels from the O stars to the WIM. Photoion- 



ization models (Hoopes and Walterbos, 2003 Wood and 



the incident photon flux. Using this fact, [Reynolds ett 
( |1995| ) found that an interstellar Lyman continuum flux 



Mathis , 2004 ) show that the spectral processing of the ra- 
diation can be complex, with the radiation between the 
H I and He I ionization edges hardening with distance 
from the source, while the spectrum at higher energies 
softens. Moreover, hot evolved low mass stars (white 
dwarfs) and interface radiation associated with the hot 
(10^~^ K) gas add harder photons to the mix (^VI). In- 
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dependently of the spectrum, the low ionization state of 
the WIM also could be the result of its low ionization 
parameter (Mathis, 1986|). 

Constraints on the fluxes of higher energy (i.e., helium- 
ionizing) photons are from observations of the He I re- 
combination line at A5876 and the [O ill] A5007 collision- 
ally excited line. Both of these transitions are promi- 
nent in O star H ll regions, where Iu{> 24 eV) is high 
enough (and is known to be high enough) to maintain 
He+ (24.6 eV) and 0++ (35.1 eV) at appreciable lev- 
els. Even qualitatively, from the flrst attem pts to de- 
tect these l ines i n the WIM ( Rey nolds[ |l985a[ [Reynolds] 
and Tufte{ [1995 ), it was clear that these ions were not 



as abundant in the WIM. More recent WHAM observa- 
tions found (He l/Ha)wiM ^ 0-5 x (He l/Ha)H ii, which 
when combined with the fact that H+/H is near unity 
(see above), implies that He+/He < 60%. The [O ill] /Ha 
results are more varied, although the ratios are typically 
less (^ 10%) those seen in H ll regions (Madsen, 200^ 



Madsen et al. 2006). The abundance of [O ill] in other 



galaxies and in the interior regions of our Galaxy can be 
signiflcantly higher than what is observed in the WIM 



near the sun dMadsen and Reynolds| |2005| |Rand[ |1997| ) . 



C. Temperature 

The temperature of photoionized gas is set by a balance 
between heating and cooling. Heat is injected by ther- 
malization of the excess kinetic energy of the electron 
duri ng the photoioni zation-recombination process (see. 

Other potential sources of heat 



Osterbrock 



1989^. 



e.g., 

could also be important, particularly at the low densi 
ties characteristic of the WIM (e.g., see Reynolds and 



Cox[ |1992j . Cooling occurs primarily from the collisional 



excitation and subsequent radiative decay of metastable 
states (i.e., forbidden lines) of the trace ions (see, e.g., 
[Osterbrock^ 1989, for a detailed discussion). The detec- 
tion of some of these "cooling lines" in combination with 
the H-recombination emission have been used to explore 
the temperature of the g as, as discussed in more detail 
below. The observations have established that 1) on av- 
erage the WIM is about 2000 K warmer than the denser, 
classical H ll regions and 2) there are signiflcant varia- 
tions in temperature within the WIM, most notably an 
increase in temperature with increasing distance away 
from the midplane, and more generally, with decreas- 
ing emission measure (or gas density). The reason for 
this temperature behavior of the WIM is not yet clear; 
it could indicate that photoionizat ion is not the only im- 
portant source of heat in the WIM (Reynolds et al. 1999 ) 



or that perhaps the spectrum of the ionizing radiation is 
modifled as it propagates far from its source (see 

Although they vary in accuracy and difficulty, three 
tools are available to explore the temperature of the WIM 
through optical emission lines: 

1. [N ll]/Ha and [O ll]/Ha trace variations in 

Tel N+, O^, and H+ are the dominant states of 



ionization for these elements in the WIM. In addi- 
tion, their emission lines have different dependences 
on temperature, so that changes in [N ll]/Ha and 
[O ll]/Ha closely track changes in Tg. Empirically, 
where < 1 R, the brightnesses of the primary 
optical forbidden lines of [N ll] , [S ll] , and [O ll] be- 
come comparable to Ha, making these lines easier 
to detect and thus attractive diagnostic tools for 
exploring variations in Tg. Calculating absolute 
temperatures is more uncertain due to necessary 
assumptions about the exact ionic fractions and el- 
emental abundances. 

2. [N ll] A5755/[N ll] A6583 measures Te directly: 

The ratio of the "auroral" (A5755) emission line 
(resulting from excitations to a metastable state 
4.0 eV above ground) to that of the much brighter 
"nebular" (A6583) transition (1.9 eV above ground) 
allows a derivation of Tg with the fewest assump- 
tions. Since this ratio involves the same ion, it is 
proportional to gC^^/^^e)^ where is the differ- 
ence in the excitation energy of the two states. Near 
8000 K, a 2000 K change in temperature produces 
about a factor of two change in the ratio. However, 
because /srss/iessa ~ 0.01, this observation is ex- 
tremely difficult for the WIM, where /hq; ~ 1 R. 

3. Widths of resolved line profiles are propor- 
tional to Ti'. Comparing the width of the Ha line 
to the width of an emission line from the heavier 
N+ or S+ ion can be used to separate the thermal 
from the nonthermal motions in the gas. However, 
very high signal-to-noise ratio line profile measure- 
ments are needed because the derived value of the 
ion temperature is proportional to the difference 
of the squares of the line widths. 

Results using these techniques in recent observations 
are summarized below. 



1. [N ii], [S ii] and [0 ii] with respect to Ha 

Two robust statements can be made about the line 
ratio observations: 

(a) [N ll]/Ha, [S ll]/Ha, and [O ll]/Ha increase with 
decreasing /hq;. In the Galaxy, this rise is most 
dramatic below /hq; = 1 R. The clearest examples 
of this are the large increases in the forbidden line 
intensities relative to Ha with increasing distance 
from the midplane, both in our Galaxy (Haffner 



eToqlBoof 



et al 1999) and others (Rand 1998 Tiillmann 



(b) Values for [S ll]/Ha and [N ll]/Ha vary greatly, but 
are strongly correlated, often with a nearly con- 
stant [S n]/[N ii] ratio over large regions. For the 
Galaxy, the ratio of [S Ii]/[N ll] does not vary by 
more than about a factor of two, except in the vicin- 
ity of a discrete ionizing source (e.g., an O star). 
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As pointed out inlHaffner et aL\ (|1999|), this behavior 



suggests that changes in the hue ratios are due primarily 
to changes in Tg. Thus statement (b) above follows from 
the fact that the [S ll] and [N ll] lines have nearly identical 
excitation energies, so that 



^6716 
^6583 



= 4.69 



(*)( 



p0.04/T4 



-0.119 



(1) 



which is only a very weak function of T4 (= Tg/lO^ K). 
From T4 = 0.5 to 1.0 with ah else constant, [S ll]/[N ll] 
decreases only about 11%. This relationship also in- 
dicates that the relatively small but real variations of 
[S n]/[N 11] that are observed in the WIM are tracing vari- 
ations of S^/N+. Combined with the very different en- 
ergies required for S+ ^ S++ (23.3 eV) and N+ ^ N++ 
(29.6 eV), we conclude that S^/S, and especially N+/N, 
vary little in the WIM and that the smaller (factor of 
two) variations in [S Ii]/[N ll] are due primarily to vari- 
ations in S^ /S. This is supported by photoionization 
models (e.g., [Sembach et al] 2000), which have shown 
that N+/N ^ 0.8 over a wide range of input spectra and 
ionization parameters. 

On the other hand, the strong temperature dependence 
of the forbidden line intensities relative to Ha is illus- 
trated by the relationship for the [N ll]/Ha 



[O 11] observations confirm that the line ratio variations 
are dominated by variations in Tg. 



2. [N II] A5755/[N 11] A6583 

One of the most direct ways of measuring Tg in ionized 
gas is to observe the ratio of two emission lines from the 
same ion but with very different excitation energies above 
ground. The /4363/^5007 ratio of [O ill] in bright H ll 
regions is perhaps the most famous of these pairs. In the 
WIM, because the [O ill] /Ha ratios are typically no more 
than 10% of that in H ll regions, the isoelectronically 
similar [N ll] line ratio, 1^755 /lebss^ has been used instead. 

By d etect ing this extrem ely weak line, [Reynolds et al 
( 2001b| ) and |Madsen| ( 2004] ) confirmed in select directions 
that Te in the WIM is indeed higher by about 2000 K 
than in the bright H ll regions. However, the details 
of the results reveal a more complicated temperature 
structure — perhaps not surprisingly. Although current 
measurements of the A5755 line still have large uncer- 
tainties. Fig. [3] indicates that Tg as inferred by the ratio 
of the [N 11] lines is systematically higher than that in- 
ferred by [N 11] /Ha in the same directions. This could 
be explained by temperature variations along the line of 
sight, since the A5755 line (excitation energy 4 eV) would 
be produced preferentially in regions with higher Tg com- 
pared to the red line (2 eV). 



H 



iV+ 
N 



no.426 -2.i8/r4 3. Line widths 



(2) 

Because N+/N and H+/H vary little within the WIM, 
variations in [N ll]/Ha essentially trace variations in Tg. 
Similar relationships can be written for [O ll], [S ll], and 
for other collisionally excited lines (see Otte et al. 2001 ). 

Using Eq. 1 and 2, we can constr uct diagnosti c dia- 
grams as pre sented by Haffner et al] ( jl999J and Madsen| 



If the intrinsic widths of emission lines can be measured 
accurately in ions having significantly different masses, 
then one can decompose the thermal motion (i.e., T^) and 
non-thermal motion contributing to their widths. H and 
S are particularly good to use because they differ in mass 
by a factor of 32, resulting in a measurable difference 
in their widths. This method has been used with some 



et al. (2006) to estimate both Tg and S^/S from obser- 
vations of l l]/Ha and [S l l]/Ha, as shown in Fig. |2] 



(adapted from Madsen 2004). With sufficient velocity 



success in both the WI M ( |Reynolds[ |1985b[ ) and H ll 
regions ( [Reynolds 1988). 



resolution it has even been possible to study variations 
in these parameters between different radial velocity com- 
ponents along the same line of sight. These results reveal 
that within the WIM there are variations in temperatures 
ranging between about 7000 K and 10 000 K and vari- 
ations in S+/S between 0.3 to 0.7. For comparison, the 
bright classical H ll regions all cluster near the lower left 
corner of the plot, [S ll]/Ha ^ 0.1 and [N ll]/Ha ^ 0.25, 
where Tg = 6000-7000 K and S+/S ^ 0.25. 

[O 11] at A3727 has a larger excitation energy than 
[N 11], making it even more sensitive to variations in 
Tg. Although inaccessible with WHAM, several extra- 



galactic studies ( jOtte et al.\ |2002[ |2001[ |Tullmann and 
[Dettmar , 2000) have traced this line, and new instru- 
ment at ionj^staxtm^to of [O 11] from the 
WIM ( .Mierkiewicz et al.\ |2006D of the Milky Way. The 



The potential power of this technique is illustrated 
in Fig. |4j which shows line width data for the large, 
high Galactic latitude H ll region surrounding the O star 
(09.5V) C Oph. The Ha from this region spans an order 
of magnitude in intensity and at the fainter end becomes 
comparable to WIM emission at low Galactic latitudes 
(~ 10 R). The ( Oph H ll region is particularly good 
for line width studies because each emission line profile 
is very well described by a single Gaussian component. 
Baker et al. (2004) measured the widths of the Ha, [S ll], 
and [N ll] emission lines from this region and derived ac- 
curate values for the temperature T^ and the mode of the 
nonthermal speeds vnt as shown in Fig. [4] The large 
ranges in these parameters appears to be real, with a no- 
ticeable gradient of increasing T^ (from 6000 to 8000 K) 
and decreasing vnt (from 8 to 4 km s~^) from the center 
to edge of the H ll region. The former is ascribed to the 
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'h« [R] [Nil] / Ha 

FIG. 2 Diagnostic line ratio diagrams. A large portion of the Galaxy in the direction of the Perseus arm (i — 130° to 160° 
and b = —30° to +30°, approximately) has been surveyed in Ha, [S ii], and [N ii] with WHAM. These panels show {left) the 
relationship between the ratios versus the intensity of Ha as well as [right) the relationship between the two ratios for the 
"local" gas component (It'LSRi < 15 km s~^). A few specific spatial regions are highlighted with grayscale fill to show effects of 
local ionizing sources: the planetary nebula (PN) S216, the H ii region surrounding the B0.5+sdO system Per, and regions 
near 0-star H ii regions. Other data points sample the WIM. In the diagram on the right, the vertical dashed lines represent 
Te, 5000 K, 6000 K, 7000 K, 8000 K, 9000 K, and 10 000 K, left to right. The slanted solid lines represent S+/S, 0.25, 0.50, 
0.75, and 1.00, lowest to highest slope. These WIM data reveal significant variations in Te and S^/S from one line of sight to 
the next. In contrast, classical H ii regions all cluster in the lower left corner of this diagram near [N ii]/Ha ^ 0.25, [S ii]/Hq; 
j^O.l. From IMadsenl [20041 



and Draine 2001 ) that begin to dominate over photoion- 
ization heating at low densities (~ 10 ~^ cm~^). 



D. Warm ionized and neutral gas 

The fact that hydr ogen i s nearly fully ionized within 
the Ha emitting gas (^II.B) implies that H+ and H^ are 

With the ad- 



hardening of th e radiation field wi th increasing distance 
from the source (Wood et al. [2005 , also, ^below), while 
the latter could be explained by a slow (2 km s~^) expan- 
sion of the H II region. A future goal is to extend this 
method to the fainter and more kinematically complex 
WIM. 

In summary, the evidence is now overwhelming from 
a variety of methods that in the WIM the temperature 
is elevated compared to the bright, classical H ll regions. 
In addition, the large variations in optical forbidden line 
strength relative to Ha within the WIM is dominated by 
changes in temperature rather than changes in ionic frac- 
tions or elemental abundances. When this is combined 
with the well-established result that the forbidden line 
ratios relative to Ha increase with decreasing Ha inten- 
sity (statement (a), above), we are led to the conclusion 
that the temperature increases with decreasing emission 

measure ( Haffner et al.^ ^1999) and thus decreasing gas |1977|)? Do the ionized and neutra l phases cycle from one 
density. These temperature variations do not appear to 
be explained solely by photoionization heating of the gas 
([Reynolds et al\ |1999[ [Wood and Mathis| |2004D, suggest- 



primarily confined to separate regions, 
vent of velocity resolved Ha surveys (i.e., WHAM), it is 
now possible to begin to explore the relationship between 
the diffuse ionized gas and neutral gas in the interstel- 
lar medium. Is the WIM the ionized portion of a low 
density "intercloud medium" (Miller and Cox, 1993) or 



is it mostly confined to the surfaces of neutral clouds, 
the transition region between cooler gas and a much hot- 



ter "coronal" temperature medium (McKee and Ostriker 



to the other? Cox and Helenius (2003) and Lockman 



(2004), for example, have suggested that portions of the 



ing additional heating sources for the WIM ([Minter and 
[Spangler , 1997| [Reynolds and Cox , 1992| Weingartner 



ionized medium may condense into neutral clouds. How- 
ever, to date there have been no observational studies of 
the H+ - H^ connection. In a general qualitative sense. 
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FIG. 3 Elevated temperature in the WIM. Select directions 
toward (□) brighter diffuse ionized regions show elevated line 
ratios in both [N ll]/Ha and [N ii] A5755/A6584 compared 
to (O) H II regions. Dashed lines mark select temperatures 
derived from each line ratio while the diagonal line traces 
"unity" in this derived temperature space. WIM directions 
have higher derived temperatures using either line ratio, but 
are also systematically above the line of "unity", suggest- 
ing that measurable temperature inhomogeneities exist along 
these lines of sight. From |Madsen[ |2004| 
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FIG. 4 Temperature Ti and the mode vnt (most probable 
value) of the non-thermal speed distribution of the diffuse H il 
region surrounding Oph. The primary axes show the mea- 
sured widths and errors of Ha and [S ll] from 126 one-degree 
pointings obtained by WHAM within a 12-degree diameter re- 
gion around C Oph. A grid of temperatures and non-thermal 
velocities has been superimposed to show the distribution of 
Ti and vnt within the H ii region. The apparent scatter in 
these values is spatially correlated (see text). The line at the 
upper-left is the line of equal widths. 



a comparison of the opti cal line profi les with the radio 
21 cm (H l) profiles (e.g., Hartmann and Burton^ ,1997| , 
indicates that at high latitudes the ionized gas tends to be 
correlated in space and velocity with the so-called warm 
neutral medium (WNM), the wide-spread, T ~ 10^ K 
phase of the H I generally associated with broad 21 cm 
profiles. There is very little correspondence between the 
Ha and the narrow-line 21 cm emission components as- 
sociated with the colder (T ~ 10^ K), denser H I clouds. 
In regions of the sky that contain anomalous velocity 
structures, specifically, the intermediate and high veloc- 
ity clouds that are not co-rotat ing with the Galactic disk, 
the correlation is quite strong (Haffner et aL\ 20011 Tufte 



et al 1998). The optical and radio emission components 



are centered at nearly the same velocities (within roughly 
5 km s~^) and have comparable velocity extents in cases 
of complicated, blended profiles. Few regions seem to 
contain only a WNM or only a WIM component. This 
relationship has been hinted at through various detailed 
absorpt ion line and emission line studies over the last 
decade (Howk et ar[|2003 Reynolds et al. |1995 Spitzer 



and Fitzpatrick[ |1993[), which together with observations 
of the [O ij A6300/Ha line intensity ratio (^II.C above) 
suggest that a significant amount of the H+ is associated 



with nearly fully ionized regions in contact with (or at 
least adjacent to) regions of warm primarily neutral hy- 
drogen. Hopefully a systematic examination of the kine- 
matic and spatial correlation between these ionized and 
neutral phases of the medium will be carried out soon. 

Although there is a close correspondence on the sky 
and in velocity between the warm neutral and ionized 
emission lines, their intensities do not appear to be cor- 
related. This has been examined in detail only toward 
two intermediate velocity H I clouds . Complexes L and K 
( |Haffner[ |2005| [Haffner et al] |200lD . In both cases, the 
column density of the neutral hydrogen Nn i and the Ha 
intensity, lua^ are uncorrelated. Whether this holds for 
more local gas in the Galactic disk still needs to be fully 
explored. A straightforward explanation for this lack of 
correlation is the fact that the intensity of the Ha is de- 
termined solely by the fiu x of ionizing rad i ation incident 
on the warm H I cloud (Reynolds et a/., 1995), which 



of course, is independent of the cloud's column density 
because the H I clouds are optically thick to the Lyman 
continuum photons. 



9 



E. The role of superbubbles 

One of the basic questions concerning the nature of 
the WIM is how ionizing photons from the O stars can 
travel hundreds of parsecs through the disk and into the 
halo. A fractal morphology of the interstellar medium 
is one possibility (see ^V|). Another is the existence of 
enormous, H l-free bubbles surrounding some of the O 
stars, which allow the Lyman continuum photons to 
travel through the cavity to ionize its distant walls (e.g., 
McClure-Griffiths et a/. ||2QQ6l [Norman and Ikeuchil [19891 



Pidopryhora et a/. | |2QQ7[ [Reynolds and Ogden[ |1979D . A 
WHAM study of one of these bubbles, the Perseus su per- 
bubble ( [Madsen €^"0/] |2QQ6[ [Reynolds ef 



2001a), has 



shown that a luminous 0-star cluster near the midplane 
can indeed produce wide-spread, nearly WIM-like ioniza- 
tion conditions out to distances of 1000 pc or more from 
the ionizing stars. However, the [N ll]/Ha and [S ll]/Ha 
ratios of the superbubble wall are not quite as large as the 
ratios observed in the surrounding WIM, suggesting that 
bubble size, gas density within the shell, supplemental 
heating, and/or the flux and spectrum of the radiation 
escaping the 0-star cluster may also be important in set- 
ting the conditions of the ionized gas. 
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FIG. 5 The dependence of DIG detections in halos of galaxies 
on the normalized star formation rate (per unit area). The 
ratio of the FIR fluxes at 60/im and 100/im (SmlS\m) is 
given versus the FIR luminosity normalized to the optical 
galaxy diameter at 25^^ mag in units of 10^*^ erg s~^ kpc~^. 
Filled symbols indicate detections of H^ gas in the halo while 
open symbols represent non-detections. Galaxies from the 
(■) starburst sample of Lehnert and Heckman[ [1995 are in- 



cluded with more normal spirals. Integrated value s for the 
Galaxy are ta ken from Bloemen et al. (1990) and Goxand 
Mezger| ( |T989| for the FIR. A Gal actlcdiskradius of 2 2.5 kpc 
is assumed. Adopted from jRossa and Dettmar] |2003a[ 



III. VIEWS FROM THE OUTSIDE: DIFFUSE IONIZED 
GAS AND STAR FORMATION RATES IN EDGE-ON 
GALAXIES 

In galaxies other than our own, the wide-spread H+ is 
most often referred to as the diffuse ionized gas, or DIG. 
While its physical properties can be measured in much 
more detail for the interstellar medium of our Galaxy, the 
detection of the DIG in other galaxie s provides the much 
needed "outside" perspective. Rand (1997, 1998^ , for ex- 
ample, was able to detect the DIG far into the halo of the 
edge-on spiral galaxy NGC 891 and found that not only 
was [S 11] /Ha anomalously high, as observed in the so- 
lar neighborhood, but that [S ll]/Ha, [N ll]/Ha, and even 
[O III] /Ha increased significantly with increasing distance 
from the midplane. This opened up a broader discussion 
about the ionization and heating processes within the gas 
and their variation with location within a galaxy. 

Furthermore, observations of other galaxies provide 
important new information about the links between the 
DIG and global properties of galaxies, such as their star 
formation activity. In particular, the structure of this gas 
perpendicular to the galactic plane, that is, in the main 
direction of the gravitational potential, is an excellent 
tracer of the dynamics of the interstellar medium driven 
by energetic galactic processes. The observations corrob- 
orate the picture of a dynamic interstellar medium driven 
by multiple and clustered supernova, producing the so- 
called disk-halo interaction between the star formation 
regions near the midplane and energized interstellar mat- 
ter that extends a kiloparsec or more above the disk. 
Early discussions of these ideas can be found in Super- 
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FIG. 6 The X-ray halo NGC 5775 from integrated EPIC 
XMM-Newton observations exhibits substructure correlating 
on galactic scales with filaments of diffuse ionized gas and 
spurs in the radiocontinuum. Adopted from |Tullmann et al.\ 
120061 



the Interstellar Medium, and the Gaseous Halo: 



The Swashbuckler^ s Approach by Heiles (1986), The disk- 
halo interaction - Superbubble s and the structure of the 
interstellar medium by Norman and Ikeuchi (1989), and 



Galactic worms by^Koo et al. (1992). 
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A. Diffuse ionized gas in the halo and star formation in the 
disk 

From an observational point of view, the study of the 
warm H+ in galactic halos — the extraplanar diffuse ion- 
ized gas — is a good start toward understanding the na- 
ture of gaseous halos and the disk-halo connection. Of 
all tracers of halo gas, the warm ionized gas is the eas- 
i est to observe with regard to sensitivity and resolution 
(iDettmarl 119921 119981 iRandl 119971 119981) , allowing us to 



study the global influence of the energy released by young 
and massive stars into the interstellar medium. 

A halo component of the DIG for external galaxies 



was flrst discovered in NGC891 (Dettmar, 1990 Rand 



et al 



galaxies (e.g 



1990| ), followed by studies of a number of other 
Pildis et 



1994 Rand, 1996 Rossa and 



Dettmar 2000 ) . In Rossa and Dettmar ( |2003a|b ) a large 
sample of edge-on spiral galaxies clearly demonstrates 
a correlation of DIG in the halo with the star formation 



rate (SFR) in the disk (see also Miller and Veilleux 2003 
for comparison). This survey covered a broad range in 
SFR extending the observations to less active galaxies. 
Until then, emphasis had been given mainl y to galaxies 
with high SFR or even star burst galaxies (Lehnert and 
Heckmai i] 1995). This relationship between the presence 
of DIG and far infrared (FIR) emission (a measure of the 
SFR) is shown in Fig. p] If the star formation rate per 
unit area is low, as determined by the FIR luminosity 
normalized to the disk surface area, the presence of halo 
DIG indeed is observed to diminish. 

The FIR luminosity per unit area therefore seems to 
be a promising indicator for the presence of halo 11+ and 
suggests that there is a minimum SFR per area required 



to drive the disk-halo interaction (Rossa and Dettmar 



2003a). Given the known large uncertainties in the nor- 
malization of supernova rates from FIR-fluxes, a reason- 
able estimate for the local break-out condition is on the 
order of ~15 SNe over the lifetime of the O star associa- 
tion. An obvious shortcoming of this analysis is the use 
of global measurements for the FIR luminosities. Future 
observations with higher angular resolution, for exam- 
ple, with the Spitzer-satellite, could allow a much better 
comparison of the local SFR with associated halo gas 
properties. 



B. The disk-halo connection and hot gas 

The scenario of the disk-halo connection predicts out- 
flows of hot (10^ K), supernova-created gas away from 
the disk and into the halo via superbubbles (or "chim- 
neys"). In this way, transparent (i.e., H I- free) pathways 
are then also provided for the transport of ionizing radi- 
ation into the halo and beyond. For galaxies with strong 
starburst activity this is well supported by observations 

oTj |2004a|b| 



(see, e.g., Cecil et al] 2002 Strickland 



et 



Veilleux et aL\ 2005| ) and described well by theory an d 
numerical modeling (e.g., Mac Low and Ferrara 1999). 
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FIG. 7 VLA radiocontinuum map of NGC 5775 at 6cm. The 
polarized intensity is given in contours on a grey scale repre- 
sentation of the Ka emission. The bars represent the magnetic 
field direction and strength. Adopted from [Tiillmann et al.\ 
|2000t courtesy M. Soida. 



However, for more normal (i.e., non-starburst) galaxies, 
we have little direct eviden ce for such outflows. High res- 
olution ground-based ( e.g., |Howk an d Savage[ |2000 ) and 
HST imaging (e.g., Rossa et a/.[|2004) studies of the DIG 
distribution in edge-on galaxies does not reveal the num- 
ber and speciflc morphology of superbubble - or chimney- 
like structures predicted by some models (jNorman and 
Ikeuchi 1989). Such structures have been found, but 
their number is small. The difficulty in flnding them 
could be an observational problem, e.g., confusion due to 



superpositions along the line of sight, or it could mean 
the the current models of the disk-halo connection need 
further work. Thus, the idea of actual mass exchange 
between disk and halo is still a working hypothesis only. 

On the other hand, there is circumstantial evidence 
for hot gas outflows in more normal disk galaxies that 
is indicated by the detection of diffuse interstellar X-ray 
emission, bo th in our Galaxy (e. g., , Smith et al. , 2007) and 
others (e.g., Wang et al. , 2001[ ) . In a sample by Ranalli 
et al. |2003") that included disk galaxies with more nor- 



mal star formation rates, a correlation was found between 
the galaxy's global X-ray emission, presumably associ- 
ated with the very hot (10^ K) gas located inside super- 
nova remnants and superbubbles, and the amount of star 
formation activity in the disk. With the current gener- 
ation of X-ray satellites it is now possible to extend the 
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X-ray studies from the starbursting and luminous galax- 
ies to less active galaxies. One result of a small survey 
with XMM-Newton is shown in Fig. je] The EPIC X-ray 
image of the galaxy NGC 5775 is integrated over all ener- 
gies and summed over all detectors. It demonstrates the 
presence of smaller scale X-ray emission features associ- 
ated with Ha filaments far above the disk. 

The diffuse H+ in galactic halos is, in most of the ob- 
served cases, also associated with synchrotron radiocon- 
tinuum halos, suggesting that there is also a relationship 
between DIG and the presence of cosmic rays and mag- 
netic fields in the halo. For example, NGC 5775 shows 
a prominent radio continuum halo, and a study of its 
polarization characteristics with the Very Large Array 
(VLA) radio synthesis telescope reveals a highly struc- 
tured large-scale magnetic field with an ordered com- 
ponent to the field lines which opens up into the halo 
(Fig. [t]). This particular structure of the magnetic field 
could be of interest for the interpretation of global dy- 
na mo theory in galax ies and is discussed in more detail 
by Beck et al (1996). The condition that allows cosmic 



rays to break out of the disk as a function of star for- 
mation rate per unit area is discussed by |Dahlem et al 
( |1995( ). A detailed study o f the kinemat i cs of t he DIG in 



NGC 5775 is presented by Heald et al (2006). A com- 



pilation of observations of edge-on galaxies in Ha, radio 
continuum, and X-rays is giv en in Dettmar (1998) and 
Rossa and Dettmar ( 2QQ3a|b ). 



In summary, studies of edge-on galaxies have shown 
that warm ionized gas halos are found in galaxies where 
the SFR per unit area is sufficiently high. Typically these 
layers of extraplanar DIG can be traced out to distances 
of z ^ 1000-2000 pc, sometimes even up to 5000 pc or 
more from the mid-plane of the disk. The H+ halos 
seem to be associated with halos of cosmic rays and X- 
ray emitting plasma, which is expected if superbubble or 
chimneys provide absorption free pathways for hydrogen 
ionizing photons as well as being conduits for the vertical 
transport of hot gas and cosmic rays away from the disk 
and into the halo. 
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FIG. 8 Deprojected profiles of the (solid) total Ha surface 
brightness into (dashed- dotted) H 11 region and (dashed) DIG 
emission for NGC 7793. From [Ferguson et a/.[[l996t 



on spirals was first noted by Monnet (1971 ), who derived 



a temperature of 7000 K, an emission measure of about 
35 cm~^ pc, and a density near 0.5 cm~^ for the emit- 
ting gas. Modern detector technology (i.e., CCDs) has 
pushed the detection of diffuse H+ to fainter regions and 
has allowed the study of other emission lines, which has 
provided insight into the relationship between the diffuse 
ionization and the O stars. 



A. Radiation from O stars and the surface brightness of 
the DIG 



IV. VIEWS FROM THE OUTSIDE: THE SOURCE OF 
THE DIFFUSE IONIZED GAS IN FACE-ON GALAXIES 

Observations of edge-on and face-on galaxies clearly 
provide complementary perspectives of the distribution 
of the DIG. Specifically, face-on galaxies show the sur- 
face brightness of emission lines from the DIG across the 
face of the galaxy, which makes it possible to explore re- 
lationships (if any) between the properties of this plasma 
and the locations and ionizing fluxes of the hot, massive 
O stars, which are the most powerful ionizing agents in 
disk galaxies and are the presumed ionizing source for 
the DIG. In general, the observations reveal a strong re- 
lationship, with the Ha flux from the DIG comparable to 
that from the classical O star H ll regions in the galaxy. 

The presence of diffuse interstellar Ha emission in face- 



Although it was not understood how Lyman contin- 
uum photons could have free paths of hundreds of par- 
sees and more in galaxies with typical interstellar H I 
densities of ~1 cm~^, ever since the discovery of diffuse 
H+, O stars have been considered the prime candidate 
for the ionization. Other known energy sources simply 
fall short in total power (e.g., Reynolds] 1984). A key 
observational step that connected the diffuse ionized gas 
to radiation from O stars was carried out by Ferguson] 
et al. (1996), who showed a quantitative relationship be- 



tween the DIG and the surface brightness distribution of 
the bright, O star H ll regions across a galactic disk. An 
example of their work is given in Fig. [8] (their Fig. 4a), 
where it is clear that the mean radial surface brightness 
profile of the DIG in Ha emission tracks that of the 
H II regions. Their study included a careful, quantita- 
tive comparison of the energetic requirements for the Ha 
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emission in the DIG with the ionizing radiation and the 
mechanical energy inferred to be emitted by the O stars 
in the H ll regions of the two galaxies measured. They 
found that the mechanical energy clearly fell short, by 
more than a factor of three. They also found that the 
maximum contribution of local sources of ionizing radia- 
tion in the DIG from stars cooler than spectral type 08 
also fell short of the luminosity required, but by a smaller 
factor. They concluded that the O star populations of the 
clusters producing the H ll regions are the likeliest main 
source of the radiation that ionizes the DIG. 

In two papers dedicated to testing the hypothesis that 
escaping Lyman continuum photons from the classical 
H II regions surrounding O stars can be sufficient to ion- 
ize the DIG, [Zurita et al\ ( |2QQ2[ |2QQQp took another im- 
portant step forward by identifying and classifying the 
H II regions in a set of photometric maps of disk galaxies 
in Ha. This process is illustrated in Fig. [9j taken from 
Zurita et al (2000), which shows in the upper left, the 



original continuum-subtracted Ha image of the nearly 
face-on spiral galaxy NGC 157, followed clockwise by a 
schematic representation of the positions and luminosi- 
ties of the classified H ll regions, a surface brightness map 
of the DIG in Ha, with the H ll regions subtracted off, 
and finally a map used for quantifying the DIG. In the 
case of this last frame, the DIG is measured by integrat- 
ing the Ha surface brightness over the full disk, using 
the values outside the H ll region boundaries and local 
mean values inside each H ll region. Some key points in 
the method for deriving these maps and quantifying the 
DIG emission include the following. First a catalog of 
H II regions was prepared, using a semiautomatic, but 
interactive method to measure their Ha luminosities, ef- 



fective radii and central positions (Rozas et a/., 1999), 



down to a lower limiting luminosity. These H ll regions 
were delimited and their emission subtracted from the 
total image by masking those pixels occupied by the H ll 
regions, but with a final refinement that contiguous pix- 
els with surface brightness higher than a set limit are also 
subtracted off. 

This procedure is illustrated in Fig. 10 fZurita et al. 



2000), which shows how this criterion for separating H ll 



region from DIG emission coincides very well with an al- 
ternative criterion in which the boundaries of an H ll re- 
gion are defined by a limiting value of surface brightness 
gradient. Having separated the H ll regions, one can then 
define the total DIG luminosity in one of three ways: (a) 
integrating the remaining surface luminosity after apply- 
ing the H II region mask; (b), as in (a) but then adding a 
contribution from the areas of the H ll regions, assuming 
this is proportional to their projected areas times their 
local DIG surface brightness, or (c), as in (a) but making 
the contribution proportional to the area of the H ll re- 
gions times the mean DIG surface brightness outside the 
regions across the disk of the galaxy. Modes (a), (b), and 
(c) give respectively lower limits, upper limits, and ap- 
proximate estimates for the total DIG emission from the 
observed galaxies. In Fig. [TT] we show the results of this 



method of estimating the DIG, for six galaxies, shown 
as the ratio of the DIG Ha luminosity to the total lumi- 
nosity for the galaxy plotted in terms of galactocentric 
radius. 

We can see that the DIG emits around a half of the 
total Ha output, that there are systematic modulations 
of this tendency with radius, and that there is slight ten- 
dency for the fraction to increase with radius. We also 
note that the projected area of the disk occupied by the 
DIG is of order 80(±10)% for all the objects shown in 
Fig. 11 In a separate study of over 100 galaxies, |Oey| 
et al. ( 2007 ) found that the amount of DIG Ha to total 



Ha from a galaxy ranged from 20% to nearly 100% with 
a mean near 60%. Voges (2006) has presented results 



suggesting an inverse correlation between the DIG Ha 
fraction and the star formation rate per unit area. 



B. An escape model for Lyman continuum propagation 

Given a full catalog of H ll region positions and lu- 
minosities for a galaxy, one can test the hypothesis that 
escaping photons from the H ll regions cause the ioniza- 
tion of the DIG by modeling the transfer of these photons 
from their points of origin. This was done in considerable 
([2QQ2| for NGC 157. This galaxy 



detail by Zurita et 



was selected because of the availability of a VLA H I 
map of reasonable resolution, as we will explain shortly. 
In Fig. [T2j we show a comparison between the observed 
surface brightness distribution in the DIG and one of the 
simplest models used. In this model 30% of emitted Ly- 
man continuum photons escape from each H ll region, 
and propagate through the DIG isotropically. 

The predicted Ha surface brightness is derived by sum- 
ming the contributions to the ionizing radiation field from 
each of the H ll regions. We can see that the result 
is remarkably similar globally to the observed distribu- 
tion, and is itself a fair verification of the initial hypothe- 
sis. However a more quantitative look at the comparison 
shows that the ratio between the predicted and observed 
DIG surface brightness is not uniform on large scales, as 
would be expected since the initial model assumes a uni- 
form slab structure for the H I involved in converting the 
Lyman continuum photons to Ha. The missing struc- 
tural parameter can be supplied by using the observed 
H I column density, as shown in Fig. [13} where we can 
see that in zones of low H I column density the ratio of 
observed to predicted DIG surface brightness is reduced. 
Maps of these two quantities give excellent coincidence 
of features, and go a step further in showing that the 
principal DIG ionization sources must be the O stars in 
luminous H ll regions. Modeling the effect of clustered 
supe r novae on the distribution of the H I, [Clarke and| 
Oey (2002| also found that the resulting dumpiness of 
the medium had a significant effect on the escape frac- 
tion of the ionizing radiation. |Zurita et al. (2002) carried 



out a number of different modeling tests of the basic hy- 
pothesis, varying the law relating the escape fraction of 
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NGC 157 




FIG. 9 Steps in the quantification of the total Ha luminosity from the DIG for the representative disk galaxy NGC 157. (a) 
Continuum subtracted Ha image. (6) Schematic form of H ii region catalog, giving position and key to the Ha luminosity 
of each region, (c) Diffuse Ha map after subtracting off the catalogued H II regions. The brightest H ii regions indicated by 
circles, (d) Measurement of upper limit for DIG. H ii regions are blanked off, then each is assigned a local value of DIG surface 
brightness. Ellipse shows limit of integrated DIG flux measured. From [Zurita et a/T]|2000| 
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NGC 157 







DIG (upper li 



EM (pc cm"^) 




FIG. 10 Technique used for separating H ii region emission 
from DIG emission, (a) Portion of Ha image of NGC 157 
with crowded field (on a spiral arm). Grey scale in emis- 
sion measure. Circles show mean catalogued radii of H ii 
regions. Dashed line at 73 pc cm~^, cut-off applied to avoid 
contamination of DIG by H ii regions. (6) Diffuse emission 
after subtracting off H ii regions and applying surface bright- 
ness cut-off. (c) Same portion of image in units of Ha sur- 
face brightness gradient. Applying uniform cut-off at 12.4 
pc cm pixel"^ (0.28" /pix) yields H II region boundaries 
equal to those found in (a), which confirms this separation 
technique, (d) Map as in (b), but with H ii region mask filled 
at level of local DIG, giving upp er limit case for tot al galaxy 
DIG luminosity (see text). From Zurita et al.\ 2000 



ionizing photons with the luminosity of an H II region, 
and varying the mean absorption coefficient of the inho- 
mogeneous neutral fraction of the DIG. However, they 
concluded that it was not possible without H I data of 
improved angular resolution to go further in testing dif- 
ferent photon escape laws, or to estimate what fraction 
of the DIG ionization could be due to mechanisms other 
than that tested. 
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FIG. 11 Radial variation of the ratio of integrated DIG lu- 
minosity in Ha to the total luminosity for six spiral galaxies 
(within O.I5R/R25 crowding precludes accurate estimates). A 
canonical value for the integrated ratio of >50% is found, with 
a total fractional area subtended by the DIG of ^-^80%. From 
[Zurita ei oT] [20001 



C. Line ratio studies 

There have been relatively few quantitative spectro- 
scopic studies of the DIG in face-on galaxies. A pioneer- 
ing study of [N 11] and Ha across the face of NGC 1068 
was carried out by Bland- Hawthorn et al (1991), who 
found very high [N Ilj/Ha ratios and discussed possi- 
ble causes for this high excitation 



Hoopes and Walter- 



bos (2003) and |Voges and Walterbos (2006) have made 
the most careful and detailed examinations to date. In 
Fig. [T4j we show the observations of the DIG close to the 
l uminous H ll region NGC 604 in the nearby spiral M33 
(Hoopes and Walterbos, 2003), obtained by placing a slit 
across the H ll region so that it sampled the DIG on ei- 
ther side. This example is representative of their study 
of the three local face-on galaxies, M33, M51, and M81. 
We can see that the line ratios, [N ll]/Ha, He l/Ha, and 
[O lll]/H/3, all tend to show higher values in the DIG. 
The full study includes measurements of these ratios as 
well as [S 11] /Ha and [S Ii]/[N ll] in H ll regions and in the 
DIG, both in the spiral arms and in the interarm zones in 
each galaxy. The results were compared with predictions 
made with "standard" H ll region ionizing fields, with 
varying ionization parameters and stellar photospheric 
temperatures between 20 000 K and 50 000 K. There 
is good agreement for the H ll regions but poor agree- 
ment for the DIG. They then considered more realistic 
radiation fields, which take into account the fact that 
the spectrum of the radiation escaping from H ll regions 
may be different from that from a pure O star (see next 
section). They modified their spectra accordingly, vary- 
ing the modeled escape fraction between 30% and 60%, 
where a lower escape fraction implies a harder spectrum. 
These results showed better agreement with the DIG ob- 
servations than the previous set, but the agreement was 



15 



NGC 1 57— Constant escape fraction (30%) 



NGC 1 57— Deprojected Hoc innage 




FIG. 12 Comparison of a DIG model with observations, (left) Modeled surface brightness in Ha of the DIG in the disk of 
NGC 157 assuming 30% of Lyc photons escape from each H ii region, and a simple propagation law through a (macroscopically) 
uniform slab model for the disk. The result is a projection in the plane of the predicted 3D Ha column density, (right) 
Deprojected image of the galaxy with H II regions masked, and a cut-off of 73 x cosi pc cm~^ applied to limit the H il 
region contamination of the DIG. This shows the basic similarity between these models and the observed DIG, though further 
refinements are important (see Fig. 13). From Zurita et al.\ 2002 



only fair. 

The authors reached the tentative conclusion that O 
star photoionization is not the sole mechanism for ioniz- 
ing and/or heating the DIG, the same conclusion reached 
by others who have studied these line ratios both in our 



Galaxy and others (e.g ., Rand 1998 Reynolds et al. 



1999 Wood and Mathis, 2004). This may indeed well be 



the case; however, before reaching a definitive conclusion, 
it would be useful t o test a modification of the types of 
models proposed in Hoopes and Walterbos (2003| based 
on the assumption that not only is the DIG itself inho- 
mogeneous, but so also are the H ll regions. The effects 
of clumping of the interstellar gas on the escape fraction 
and the spectrum of the ionizing radiation is discussed in 
more detail below. 



V. MODELING THE WIM/DIG: EFFECTS OF 
RADIATION TRANSFER THROUGH A CLUMPY 
INTERSTELLAR MEDIUM 

Models of a more realistic (i.e., 3-D clumpy) medium 
may provide a clearer view of the ionization and heating 
processes in the gas as well as insight into such questions 
as: 

• What are the constraints on the structure of the in- 
terstellar medium that allow photons to penetrate 



from the midplane O stars to large distances in the 
halo? 

• How much of the Ha comes from the ionized 
surfaces of dense clouds and how much from a 
smoother low density medium occupying the space 
between the clouds? Is a cloud-intercloud model 
sufficiently realistic? 

• How much extra, non-ionizing heating is necessary 
to explain some of the line ratios? 

• What is the role of interfaces between the ionized 
and neutral gas? Are photoionization codes missing 
crucial physics at interfaces? 

• How much ionizing radiation escapes from galaxies, 
and what is the spectrum of the escaping radiation? 



A. What exactly is an H ii region? 

Before describing various modeling efforts for the DIG, 
we will try to envision the structure of an H II region. In 
the classical picture of a Stromgren sphere, a spherical 
volume of interstellar gas of uniform density is ionized 
by a central source. For plausible values of the ioniza- 
tion parameter there is a s harp transition from H ll to H I 
at the Stromgren radius (Osterbrock, 1989 Stromgren 
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FIG. 13 Ratios of observed to modeled DIG of NGC 157, for 
three variant models of photon escape from H ii regions, as- 
suming a slab structure to the H i in the disk, (a) Constant 
fraction of photons emitted by the O stars escape from each 
H II region. (6) Only H ii regions with luminosities higher 
than a critical value show significant Lyc photon escape, (c) 
A constant underlying escape fraction for all H ii regions plus 
an increasing increment for H II regions above a set Ha lu- 
minosity. The results are similar for all three, (d) Observed 
H I column density map (with catalogued H II regions over- 
laid in black). Note complete coincidence of low zones of 
observed/modeled ratio with zones of low H i column density, 
as predicted if photon escape from the H II regions is the main 
ionizer of the DIG. From ,Zurita et aE] [20021 
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Distance Along Slit (pc) 

FIG. 14 Selected emission line ratios from long slit spectra 
across the luminous H II region NGC 604 in M33. The solid 
line is the Ha surface brightness showing rising to the center 
of the H II region, and falling to low values in the DIG outside 
it. Note the tendency of the line ratios to rise in the DIG. 
From [Hoopes and Walterbosl [20Q3| 



a more realistic picture of an H ll region comprises high 
density clumps embedded in a lower density interclump 
gas. 

Depending on the density of the clumps they may be 
fully ionized or be dense enough so that they hav e an 
ionized skin and a shielded, neutral interior (e.g 



ammanco et al 



qualitatively exp 



Gi- 



2004). This more complex picture can 



ain the observed features of H ll regions 

including ^^blis ter" regions such as the Orion Nebula (e.g., 

1939| ). However, in order to better explain line ratios in Ferland, 2001). Consequences of this clumpy structure 



real H ll regions, the "filling factor" approach was devised 
in which the ionized gas is confined to small, fully ion- 
ized clouds surrounded by vacuum, and with the ionized 
gas o ccuping only some fraction of the H II region vol- 
ume dOsterbrock and Flather , 1959 Osterbrock , 1989 



for emission line ratios from H ll regions were quantified 
in 



Giammanco et al. (2004) and in Giammanco et 



( 2QQ5| ). Based on the models explained in these two arti 



Stromgren[ 1948^. Both of these pictures deviate from 
real H ll regions, where neutral condensations can exist 
within the H ll region. In addition there may be a con- 
voluted interface between ionized and neutral gas on the 
faces of clouds exposed to ionizing radiation. These in- 
terfaces and the shadowed regions behind the clumps will 
have a very different ionization and temperature struc- 
ture and hence also a different emission line spectrum 
compared to fully ionized clouds or the gas inside a uni- 



form Stromgren sphere (e.g., Williams 1992). Therefore, 



cles. Fig. 15 shows an example of two spectra of ionizing 
photons emerging from an H II region, one hardened by 
partial ionization of H I while propagating through the 
denser clumps within the region, and the other essen- 
tially unaffected as it has passed through the interclump 
gas. Applications to the DIG will entail the combination 
of these two types of spectra with suitable weightings. 
When this has been done, we will be in a position to 
make clearer statements about what fraction of the ion- 
ization and heating of the DIG can be due to photons 
escaping from H II regions, the presumed primary source 
of the ionization. 
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FIG. 15 Components of unattenuated and hardened ionizing 
spectra emerging from a clumpy model H II region, with ge- 
ometrical filling factor of denser clumps 10~^, clump density 
100 cm~^, and mean clump radius 1 pc. Solid line: input spec- 
trum emerging diluted but unchanged in relative intensities 
through interclump medium. Dashed line: hardened spec- 
trum after transiting internal clumps. In this model the DIG 
would be ionized by a combination of both spectra, weighted 
according to the details of the clump model. For details of 
these models see |Giammanco et al.\ |2004[ 



B. Modeling the ionization structure of the DIG 

A complex 3D picture is almost certainly needed to 
explain the DIG. As we will see below, for example, a 3D 
medium can help to answer questions regarding the pen- 
etration of ionizing photons to large distances. However, 
among the most crucial parts of any 3D photoionization 
model are the density distributions of its gaseous com- 
ponents. If, for example, we adopt a model consisting 
of dense H I clouds within a uniform lower density in- 
tercloud H I medium, then we will need to know what 
is the density of the intercloud medium, how much does 
it contribute to the observed Ha, and what is the cov- 
ering factor of the clouds? In models with more com- 
plex density distributions, we would need more accurate 
knowledge about the actual spatial and density distribu- 
tions of the ionized as well as the neutral gas, which are 



2008 Kim 



only b eginning to be explored (e.g.. Hill et al. 

2008j). Questions pertaining to the density struc- 



et al. 



ture of the interstellar medium are the major unknowns 
in photoionization models. 

Until recently, multi-dimensional photoionization 
models of the DIG did not consider the faint forbid- 
den line emission, and focussed mostly on Ha and the 
3D ionization structure of hydrogen via "Stromgren vol- 



ume" techniques (Dove and Shull, 1994; Dove et al.^ 2000 


Franco et al. 


1990 


Miller and Coxl 1993) or Monte Carlo 


simulations ( 


Ciarc 


i et al. 2002 


Wood and Loebl 2000). 



These models demonstrated that extended DIG layers 
could be produced by the ionizing radiation of O stars 
in a clumpy medium, where the 3D density structure of 



the medium was simulated using a variety of mo dels in- 
cluding "standard clouds" (Miller and Cox, 1993), super- 
bubbles formed by the action of supernovae and stellar 



winds ( Dove et al. 
sities (Ciardi et al. 



2000), and two-phase and fractal den- 



2002 Wood and Loeb, 2000). A 3D 



density structure was found to be necessary to allow ion- 
izing photons to penetrate to large distances from the O 
stars, which are confined to discrete locations near the 
galactic midplane. In a smooth medium with a typi- 
cal interstellar density of 1 cm~^, an O star will form a 
Stromgren Sphere of radius about 60 pc, and within the 
neutral medium an ionizing photon will penetrate only 
0.1 pc; so, 3D structures must provide the low density 
paths that allow the photons to traverse kiloparsec size- 
scales and ionize the gas far from the O stars and at large 
heights above the midplane. 



1. Escape of ionizing radiation through superbubbles 



The work of Dove et al. (2000) investigated whether 



superbubbles created by the dynamical action of super- 
novae and stellar winds could provide the low density 
paths. They found that the dense swept up shell of ma- 
terial would trap ionizing photons within the bubble and 
thus inhibit the escape of ionizing radiation to the halo, 
unless, of course, the size of the bubble were so large 
that it reached into the halo. Their models employed a 
smooth, continuous structure for the swept up shell. The 
fragmentation of shells and the punching of holes by su- 
pernovae and stellar winds as they expand into a clumpy 
medium may offset the trapping of radiation in smooth 
models and allow a larger fra ction of ionizing phot ons 
to escape. Recent analysis by Terebey et al. (2003) of 



the W4 chimney suggests that a significant fraction of 
ionizing radiation is able to escape from a fragmented 
shell-like structure near the midplane, perhaps to be ab- 
sorbed on the distant wall of a much larger cavit y that 
reaches into the lower halo ( [Reynolds et al. , 2001a). Pho- 
toionization models of bubbles created by stellar winds 
in a turbulent star forming cloud also suggest that large 
escape fract ions (in excess of 20%) are possible ( [Dale 
et al. 2005). Future work should study in more detail 



the ionization structure of 3D hydrodynamic simulations 
of supernova created bubbles and superbubbles to deter- 
mine escape fractions and the spectrum of the escaping 
ionizing radiation. A comparison of such models with ob- 
servations of H II regions and the surrounding DIG would 
allow a critical test of the predicted interstellar structure. 



2. Two- and three-dimensional ionization structure of the DIG 



The [Miller and Cox| ( [19931 ) model of the DIG adopted 
a two component, vertically stratified density structure 
n(H) = 0.1exp(-|z|/0.3) + 0.025exp(-|z|/0.9), where 
the number density is per cm^ and the z-distance is 
in kiloparsecs. This represents the concentrated neutral 



18 



layer and the extended ionized layer. They also included 
an approximation for absorption by dense opaque clouds 
using a model that reproduces the statistics of clouds 
in the interstellar medium. Although their smooth den- 
sity is lower than the average density inferred for and 
H+, when implemented in their "standard cloud" model 
including the known ionizing sources in the solar neigh- 
bourhood, it reproduced the average observed emission 
measures (EM) and dispersion measures (DM). 

An obvious criticism of the Miller & Cox model is that 
the mean density they used is smaller than that inferred 



for th e H I in the Galaxy (e.g., [Dickey and Lockman 



1990). However, they did use the known distribution 



and ionizing luminosity of O stars in the Solar neigh- 
borhood and showed that the above density distribution 
would allow for the gas to be ionized to large \z\. In re- 
ality the gas is clumped on a wide range of size scales, 
so the Miller & Cox density provides an estimate of the 
maximum amount of smooth component allowed, if ion- 
izing photons are to penetrate to large \z\. Therefore, 
the Miller & Cox density profile may be a good starting 
point for the smooth intercloud medium in 3D models 
that incorporate both smooth and clumpy components. 
To illustrate the effect of a clumpy, 3D density distri- 



bution. Fig. [16] shows slices through the ionized volumes 
created in a completely smooth medium and in a medium 
with both clumpy and smooth components. In each case 
to simulate an O star, the ionizing source is a 40 000 K 
blackbody point source at the center of a grid which is 
2000 pc on a side. The hydrogen ionizing luminosity Q is 
varied in the range 10^^ s~^ < Q < 10^^ s~^, thus giving 
a range of sizes for the ionized volumes, corresponding to 
that for a single O star up to that for a large O star as- 
sociation. For the smooth mediu m the density structure 
is that of a Dickey-Lockman disk (Dickey and Lockmanj 
P-990J . For the clumpy medium the average |z|-height de- 
pendent density is also that of a Dickey-Lockman disk, 
except that the gas has been turned in to 3D clumps usin g 
the hierarchical clumping algorithm of Elmegreen| ( 1997 ) . 



For this simulation one third of the mass remains in a 
smooth component and the remainder is converted into 
hierarchical clumps. The photoionization simulation is 



performed using the 3D Monte Carlo code of |Wood et al 
( [2004, ). 

The contours in Fig. [16] show slices through the ion- 
ized regions for different source luminosities. It is im- 
mediately seen that for a given source strength, the 3D 
density structure with clumps allows ionizing photons 
to reach much larger \z\ distances than the completely 
smooth distribution. In the smooth model the source 
with Q = 10^^ s-^ (typical of a single O star) creates a 
very small ionized volume, while the same source ionizes 
a much larger volume in the 3D simulation. Recall that 
both the smooth and clumpy density grids have the same 
mean n{z) given by the Dickey-Lockman density distri- 
bution. Therefore, it appears that 3D density structure 
readily allows for ionizing photons to penetrate to large 
1 2; I, thus solving a problem common to models that adopt 





FIG. 16 Slices through the ionization structure for point 
sources in (left) smooth and (right) hierarchically clumped 
density distributions. The mean vertical density structure for 
each simulation is that of a Dickey-Lockman disk. From the 
inner to outer contours the source luminosities (photons per 
second) are 10^^ 3 x 10^^ 5 x 10^^ and 10^°. Comparing the 
contours in the smooth and clumpy models, it is clear that 
the 3D density structure provides low density paths allow- 
ing photons to reach much larger distances from the ionizing 
source. 



a smooth density structure of the medium. 

C. Modeling the emission line spectrum of diffuse ionized 
gas 

The physical conditions in the DIG are revealed 
through its emission line spectrum, which probes its ion- 
ization state and temperature structure (^Il|. The main 
observational characteristics that successful models of the 
DIG spectrum must address are: 

• The temperature, probed via the dominant ion ra- 
tios [N 11] /Ha and [O ll]/Ha, appears to rise with 
increasing |z|-height above the plane. Most of the 
nitrogen and oxygen are N+ and 0+, and [N ll] 
and [O 11], excited by electron collisions, are very 
sensitive to temperature. 

• Compared to traditional H ll regions the DIG ap- 
pears to be ionized by a softer radiation field. This 
is revealed in the elevated [S Ii]/[N ll], a measure of 
S+/S, and the apparent underionization of He im- 
plied by low He l/Ha recombination line intensity 
ratios. 

• WHAM observations show that [O l]/Ha is weak, 
at least near the Galactic midplane, while in 
NGC891 [O l]/Ha is observed to increase with in- 
creasing distance above the plane. Models predict 
that interfaces between the ionized and neutral re- 
gions should be revealed by increased [O l]/Ha and 
a rapid rise of [N ll]/Ha and [S ll]/Ha near the 
transition region, where the gas temperature is ris- 
ing rapidly. 

• In some galaxies [O ill] /Ha is seen to remain high 
or even increase at large |z|-heights, whereas pho- 
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toionization models predict [O ill] /Ha to decrease 
with height above the plane. 



1. One-dimensional models 

While the 2D and 3D techniques described above were 
used to explain the formation and structure of the DIG 
via 0-star photoionization, detailed ID photoionization 
codes were used to explain the DIG's emission line spec- 



trum (e.g., ^Collins and Rand 


2001 [Domgorgen and] 


Mathisl|1994| Hoopes and Walter! 


bos| 20031 Mathisl 1986| 


2000 Sembach et a/., 2000). In general, the ID sim- 



ulations parameterized the problem with an ionization 
parameter, calculated the ionization structure of a con- 
stant density spherically symmetric volume, and formed 
line ratios by taking ratios of the total emission line in- 
tensities from the entire volume. By combining simu- 
lations with different ionization parameters and differ- 
ent supplemental heating rates, these volume averaged 
models could explain many of the observed line ratios in 
the DIG. However, ID models do not allow for a self- 
consistent explanation of the observed trends of line ra- 
tios. In particular, [N ll]/Ha and [S Ilj/Ha are observed 
to be anticorrelated with the Ha intensity and increase 
with increasing distance from the midplane (e.g., Haffner 



et a/.||1999||Rand||1998[ ). 

Explaining these observational trends requires models 
that do not use volume averages, but form line ratios 
for lines of sight that pierce through the ionized vol- 
ume at increasing impact parameters away from the ion- 
izing source. The pl ane parallel slab models presented 
by Bland- Hawthorn et al.\ ( 1997) showed how the depth 



dependence of the ionization and temperature structure 
(see their Fig. 9) could explain the trends for [N ll]/Ha 
line ratios with distance fr om the source. More recently, 
Elwert and Dettmar ( |2005 ) presented a grid of plane par- 
allel models that could explain the observed increase of 
[N ll]/Ha and [S ll]/Ha above the midplane. 

A persistent problem with all these photoionization 
models is that they cannot explain the rise of [O ill] /Ha 
with height above the plane unless some other ionization 
and/or heating source is invoked, such as shoc ks, photo- 
electric heating, or turbulent mixing layers (Collins and 



Rand 2001 Reynolds et al] |1999[ SlavhieFo/T 1993 ) . 



Also, unless the medium is fully ionized beyond the edge 
of the simulation (i.e., density bounded) the models pre- 
dict uncomfortably large [O l]/Ha ratios at the edge of 
the ionized volume. Clumping in H ll regions may explain 
the apparent absence of these edge effects as discussed 
below. 
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FIG. 17 Ha and [N ii]/Ha maps for (upper panels) uniform 
density and (lower panels) hierarchically clumped H II region 
models. In each simulation a uniform background has been 
added to the Ha and [N ii] maps to simulate the H ii region be- 
ing surrounded by low density DIG. Notice the very irregular 
boundary in the clumpy model and that the rise of [N ii]/Ha 
with distance from the central source is much steeper in the 
smooth model. 



ionization and temperature structure and emission line 
strengths for multiple ionizing sources within 3D geome- 
Wood and Mathis (|2004) presented simulations 



tries. 

for photoionization and line intensity maps of a strati- 
fied interstellar medium. Their models showed that due 
to hardening of the radiation field in the spectral region 
between the H I and He I ionization edges, the temper- 
ature naturally increases away from the ionizing source. 
They also found that [N ll]/Ha and [S ll]/Ha rise with 
\z\ in their vertically stratified density structure and that 
[S Ii]/[N ii] remains fairly constant with \z\ as observed. 
While their models could explain most of the line ratios, 
additional non-ionizing heating was still required to pro- 
duce the largest line ratios observed. 



2. Two- and three-dimensional models 



3. Interfaces and three-dimentional H ii regions 



Recent advances in the development of Monte Garlo 



photoionization codes ([Ercolano et al. 2003 Och et al. 



1998 Wood et al., 2004) now enable calculations of the 



In standard photoionization models of the interface be- 
tween ionized and neutral gas, the medium is turning 
neutral and the temperature is rising rapidly due to the 



20 



hardening of the radiation field (e.g., [Wood and Mathis 



[2QQ4J . For optical observers, such an interface should be 
characterized by enhanced emission of [O l], [N ll], and 
[S ll] relative to Ha. As mentioned above these effects 
have in general not been observed in the DIG. However, 
recent results from the WHAM survey of low emission 
measure H ll regions do appear to show an increase of 
[N ll]/Ha and [S llj /Ha away from the ionizing sources 



(Wood et fl/., 2005). Fig. 17 shows that models of hi- 



erarchically clumped H ll regions give a shallower rise 
of [N ll]/Ha away from the ionizing source than that 
predicted by uniform density models. This is because 
the edge of the H II region is very irregular and sight- 
lines probe different columns of ionized gas and different 
temperatures compared to the regular, rapid rise in tem- 
perature and neutral fraction seen at the edge of a uni- 
form density H ll region. This shallower rise is in better 
agreement with the observations. Further work extend- 
ing 3D models from individual H ll regions to study pho- 
toionization in global models of the interstellar medium 
is required. This will test whether interfaces can indeed 
be masked by 3D effects and the edges of H ll regions 
may be lost in diffuse foreground and background (non- 
interface) emission, although the observed weakness of 
the [O l] emission still presents a problem. Most pho- 
toionization codes neglect the effects of shocks and the 
expansion of ionization fronts at the edges of H ll regions. 
Perhaps there is some crucial missing physics from cur- 
rent photoionization codes that is important for modeling 
interface emission. There are efforts underway to include 
such effects in the CLOUDY photoionization code (Elw- 
ert 2005, private communication). 



4. Leaky H ii regions and the He^/H^ problem 



Leaky H ll regions appear to be an important source 
of ionizing photons for the DIG (see the discussion in 
^IV). The ionizing photons escape from them either di- 



rectly through empty holes or via escape through density 
bounded ionized gas. In the first case the escaping spec- 
trum will be that of the ionizing star, while in the second 
case the spectrum will be modified due to partial absorp- 
tion of the radiation by neutral H and He. Compared to 
the intrinsic stellar spectrum the escaping spectrum is 
generally harder in the H-ionizing continuum and has its 
He-ionizing photons supressed (e.g., ^Hoopes and Walter-, 
bosl 120031 iWood and Mathisl 120041). 



A combination of direct escape and transmission 
through density bounded ionized gas is likely to occur 
in a 3D H ll region. In fact, line ratios in the DIG 
surrounding leaky H II regions may allow us to deter- 
mine whether the escaping photons are dominated by 
those escaping through true holes in the H ll region or 
by those penetrating density-bounded H II gas. Specifi- 
cally, the supression of He-ionizing photons penetrating 
density bounded H ll regions may explain the underion- 
ization of He in the DIG. However, more observations of 



He I and other lines as a function of distance from H ll 
regions in our Galaxy and others are required to test this. 
It cannot yet be ruled out that most of the DIG ionization 
is produced by later O stars (which have a soft spectrum 
consistent with the He I observations), because the earli- 
est (hard-spectrum) O stars may preferentially have their 
ionizing radiation absorbed within their parent molecu- 
lar clouds. However, a study of face-on galaxies suggests 
that th ese cooler O stars do n ot provide sufficient ionizing 
power (Ferguson et al. 1996). Clearly more observational 
and theoretical investigations are required into the role 
of leaky H II regions and the spectrum of their escaping 
ionizing radiation. 



VI. IONIZING RADIATION FROM HOT GAS-COOL GAS 
INTERFACES 

Early-type stars produce a prodigious amount of ion- 
izing radiation and are capable of ionizing gas well above 
the midplane of the Galaxy. But does this stellar radi- 
ation explain all of the observed H+? For example, do 
some of the O star photons completely escape the Galaxy 
to account for the ionized gas observed in the high veloc- 
ity H I clouds (HVCs) and the Magellanic Stream, located 



up to 50 000 pc from the Galaxy (Putman et al. 
Tufte et all 119981 iWeiner and Wi 



liamsl 1996)? If so. 



2003 



then the Ha surface brightness of these clouds provides a 
direct measurement of the flux of Lyman continuum ra- 
diation that completely escapes the Galaxy (e.g., [Tufte 
1998) . However, because these cooler clouds ap- 



et al. 



pear to be immersed in a much hotter plasma (e.g.. Sav- 



age et al. 2003), we must at least consider the fact that 



hot gas-cool gas interfaces are also a source of ionizing ra- 
diation. Such radiation may even play an important role 
in the ionization of a very local H I cloud in the vicinity 
of the sun. 

Observations of 0+^ a nd detections of other ions in 
high states of ionization (Savage et al. 
and Savage, 1992; Sembach et 
faces between hot ( 



2000 Sembach 



2003 ) show that inter- 
10^ K) g'^s and cooler (< 10"^ K) 
gas are widespread throughout the interstellar medium 
and Galactic halo. The resulting intermediate tempera- 
ture (~ 10^ K) gas at these interfaces produces extreme 
ultraviolet ionizing radiation. Thus, even though O star 
photons can leak through a clumpy interstellar medium 
and/or through superbubble chimneys, interfaces have 
the advantage that they exist wherever the hot gas and 
cooler clouds do, including places where ionizing radia- 
tion from O stars does not reach. Although significantly 
weaker than the flux from hot stars, interface radiation 
may be more widely distributed, and because the emis- 
sion is generated in a thin layer adjacent to the absorbing 
cloud, interface radiation is efficiently used for ionizing 
that cloud as well as being a source of ionization for other 
more distant clouds. 
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A. Types of interfaces 



A variety of different types of interfaces may exist, de- 
pending on the physical processes operating and dynam- 
ical state of the boundary region between the hot and 
cooler gas. These include : evaporative boundaries (e.g., 
|Cowie and McKee, 1977 ), cooling/condensation fronts 



2.0 



dShapiro and Benjamin[|1991 ), and turbulent mix i ng lay - 



ers ( [Begelman and Fabian 



19901 Slavin et a/., 1993). 



In evaporative boundaries, thermal conduction heats the 
cool cloud and produces an outflow. This requires that 
the magnetic field topology is such that the warm gas is 
not shielded too thoroughly from the hot gas. In cool- 
ing/condensation fronts, slow accretion of hot gas onto 
the cool gas occurs as the hot gas cools radiatively. A tur- 
bulent mixing layer (TML) can develop in regions where 
there is shear flow at the hot / cool boundary that leads to 
hydrodynamical instabilities and mixing of the hot and 
cool gas. The mixed gas in a TML cools rapidly due 
to its temperature and non-equilibrium ionization state. 
Although all these types of interfaces share some charac- 
teristics, the ionization state of the gas can be radically 
different for different types of interfaces. For example, 
relative to collisional ionization equilibrium the gas can 
be highly underionized in evaporative outflows, overion- 
ized in cooling/condensation fronts, or a combination of 
both overionized and underionized as in a TML, wherein 
the formerly cool gas is underionized and the formerly hot 
gas is overionized. In Fig. [18] we show the results of pre- 
liminary numerical hydrodynamic simulations of a TML. 
Other possibilities exist for hot gas/cool gas interfaces in- 
volving various combinations of cooling, conduction and 
mixing, but those have yet to be explored. 

The flux and spectrum of ionizing radiation emitted, 
is determined by the ionization-temperature-density pro- 
flle in the interface. In general, underionized gas radiates 
more strongly than overionized gas, because ions that are 
being ionized up will generally be excited several times 
before being ionized. In Fig. [19] we show a comparison of 
the EUV/soft X-ray spectrum generated in an example 
evaporating cloud boundary and TML. In this calculation 
the hydrogen ionizing photon production rate between 14 
eV and 24 eV is approximately 2 x 10^ photons cm~^ s~^. 
This is only about 10% the ionizing flux that appears to 
be in cident on high velocity clouds in the Galactic halo 
(e.g., Tufte et al] 1998, 2002), for example; however, the 



uncertainty in the properties and morphology of actual 
cloud interfaces (e.g., the number of interfaces a line of 
sight through a cloud intersects; see Fig. 18) leaves open 



the possibility that the Ra produced by interface radia- 
tion could be signiflcant in regions where stellar ionizing 
photons do not reach. 



B. A test case 

The Local Interstellar Cloud (LIC) that surrounds the 
Solar System appears to be an excellent candidate for 
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FIG. 18 Temperature in a shear layer between hot and warm 
gas. The hot gas {top, T = 10^ K) flows to the left at 10 
km s~"^ while the warm gas {bottom, T = 8000 K) flows to 
the right at 10 km s~^. Cooling is not included in this 2-D 
calculation. 
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FIG. 19 Comparison of ionizing flux generated in evaporating 
cloud with that generated in a turbulent mixing layer. 



exploring interstellar interfaces. It is inside the Local 
Bubble ( |Cox and Reynolds 1987), and thus probably 
surrounded by hot gas. There is no direct Lyman con- 
tinuum radiations from O stars and the nearest B star (e 
CMa) is more than 100 pc away. The ionization of the 
LIC is very well characterized (it is best observed inter- 
stellar cloud in the Universe), and its ionization state is 
somewhat unexpected, that is, quite different from the 
WIM. In the LIC, the hydrogen is only moderately ion- 
ized at ~20-40%, and He is more ionized than H. 
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Models that inc lude ionizing radiation fr o m an evap- 
orative interface (Frisch and Slavin, 2003 Slavin and 



Frisch 



2002), are generally successful in matching the 



myriad of data available and indicate that a diffuse EUV 
source above the weak ionizing flux provided by nearby 
stars is necessary to explain the ionization. Problems 
remain, however, in explaining the relatively low column 
densities of 0+^ and as well as the high column of 



Si+^ that have been observed (Gry and Jenkins, 2001 



Oegerle et al, 2005). These discrepancies would seem 



to point to a different type of interface surrounding the 
Lie. The existence of the LIC does raise the question of 
how common such partially ionized clouds are and how 
much they contr ibute to the diffuse interstellar H+ (e.g., 
Reynolds, 2004). While the density and temperature of 



the LIC are very close to what is found in the WIM, 
the low values of [O I] /Ha observed by WHAM imply 
that H+/H for most of the diffuse ionized gas in the so- 
lar neighborhood is much closer to unity than the lower 
value found for the LIC (although He+/He may be simi- 
lar). 

In summary, while many aspects of the physics of inter- 
faces are yet to be explored, the fluxes produced by such 
interfaces are probably weaker than the ionizing fluxes 
require to produce the WIM (^Il|). Also, the ionization 
state of the LIC cloud suggests that such clouds can ac- 
count for only a small portion of the ionization associated 
with the WIM. Nonetheless, there are several conclusions 
we can draw from existing observations and theory that 
are relevant to the WIM: 

1. Hot-cool gas interfaces appear to be widespread 
through the interstellar medium. Wherever in- 
terfaces with intermediate temperature gas (T ~ 
10^ — 10^"^ K) exist there will be ionizing radiation 
generated, producing diffuse ionization. The fact 
that the spectra peak in the EUV and that a large 
fraction of the generated flux should be captured by 
the cool cloud gas makes interface radiation a po- 
tentially significant source of ionizing radiation in 
regions shielded from direct 0-star Lyman contin- 
uum photons. For example, if the interface is very 
convoluted, a typical line of sight through a cloud 
may pass through many hot-cool surfaces, result- 
ing in a total Ha intensity for the cloud that is 
comparable to what is observed for the HVCs. In- 
terface radiation appears to be require to explain 
the properties of the local cloud near the sun. 

2. The theoretical spectra from evaporating clouds 
and turbulent mixing layers are harder than that 
of the O stars, which could result in the presence 
of some ions in the WIM that O stars cannot pro- 
duce. 

3. The type of interface that exists at the boundary of 
a cloud depends (at least) on the dynamics of the 
boundary, the geometry of the cloud and the degree 
to which thermal conduction occurs. Thus refine- 
ments in theory and observations of good candidate 



clouds are needed to determine the extent to which 
interfaces contribute to ionization in our Galaxy 
and others. 



VII. SOME QUESTIONS FOR FUTURE STUDY 

A. What is the source of the elevated temeratures? 

Some of the line ratios, including [O ill] /Ha and the 
largest observed [N ll]/Ha and [S ll]/Ha ratios do not 
seem to be explained by photoionization, even when 
hardening of the O star's radiation is considered (fv|). 
Do shocks pla y a role (Collins and Rand 2001 Hidalgo- 
Gamez , 2005[ )? Or is there some additional source of 
no n- ionizing heat, such as photoelectric heating from in- 
terstellar dust particles or large molecules (Reynolds and| 
Cox 1992 Weingartner and Draine[^2001 ), dissipation 
of turbulen ce ([Minter and Spangler " p!997 ), magnetic re- 
connection ( Raymond]" 1992 [), shocks an d cooling hot gas 



(Collins and Randl 20011 Slavin et al 1993)? Most of 



these mechanisms seem capable of generating heating 
rates of the order 10~^^ ergs~^ cm~^, the power require- 
ment for the additional heat source (see discussion in 
Reynolds et al. ,^1999). In addition, , Reynolds et al. (1999) 
show that for these non- ionizing sources, their heat can 
dominate over that from photoionization at low (< 10~^ 
cm~^) densities, because their heating rates are propor- 
tional to the first power (or less) of the density, rather 
than the second power as with photoionization. This 
could explain the observed inverse correlation between 
the line ratios and the Ha intensity (see Fig. [2|. Ways of 
discriminating between these mechanisms are very much 
needed, but none is yet forthcoming. 

Can photoelectric heating be eliminated as a candidate 
for the supplemental heat? In neutral gas, the photoe- 
jection of electrons from polycyclic aromatic hydrocarbon 
(PAH) molecules by stellar ultraviolet photons below the 
Lyman limit is believed to be a major source of heating. 
However, PAHs may not be abundant in ionized g as, as 
suggested by both observational and theoretical consider- 
ations. Observati ons and modeli ng of the Orion Nebula 
H II region (e.g., Ferland, 2001), show that PAHs are 



not present in the ionized gas there. From a theoretical 
perspective, if PAHs were present in H ll regions their 
large Lyman continuum opacity (e.g., Weingartner and 
Draine, 2001) would compete with hydrogen for ionizing 



photons and result in around 90% of the ionizing pho - 
tons being absorbed by PAHs (Mathis and Wood, 2005|. 
PAH opacity peaks around 17 eV, which, for example, 
would result in He^/H+ ratios in ionized gas containing 
PAHs being much higher than is observed in H ll regions 
surrounding late O stars and in the DIG. On the other 
hand, silicaceous dust particles have a grayer wavelength 
dependent opa city and only absorb abo ut 5% of ionizing 
photons in the Mathis and Wood (2005) simulations, not 



significantly changing the He+ abundance or other line 
ratios. However, the role they may play in photoelectric 
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heating is unclear. 



D. Do hot, pre-white dwarf stars play a role? 



B. What is the spatial distribution of the gas? 

Does the wide-spread Ha emission originate from a low 
density intercloud gas that occupies large volumes of the 
interstellar medium? Or is some significant fraction of 
the Ha coming from the denser ionized faces of clouds 
or superbubble walls? Models investigating the differ- 
ent scenarios and comparisons of the resulting emission 
line and line ratio maps with observations, including elec- 
tron column densities from pulsar dispersio n measures 
may provide some answers. Hill et al. (2008), for exam- 
ple, have shown that the statistical distribution of emis- 
sion measures in the WIM is matched very well by the 
density distribution produced in magnetohydrodynamic 
(MHD) models of an isothermal, mildly supersonic tur- 
bulent medium. Also, new observational methods that 



sample the H+ at very small scales (Hill et al. 2003) or 



that access emission lines in a different spectral region 
(e.g., Kutyrev et a/., 2004|, may offer new insights about 



the small scale structure and dynamics of the gas as well 
as its larger scale distribution within the disk and halo. 



In their late stages of evolution, low-mass stars pass 
through a hot photospheric phase after shedding their 
outer envelopes. Once the stellar envelope (i.e., the plan- 
etary nebula) has expanded sufficiently and is optically 
thin to the ionizing radiation from the stellar core, the 
Lyman continuum photons are available to ionized the 
ambient interstellar medium. Such stars may be respon- 
sible for the ionization of small, localized regions in the 
low density interstellar medium (Reynolds et al 2005), 
but it is not known whether they also contribute sig- 
nificantly to the more diffuse ionizing radiation produc- 
ing the WIM. Their luminosities and their lifetimes in 
this phase are orders of magnitude smaller than that 
of massive O stars; however, in comparison to O stars, 
their numbers are enormous, and they are much more 
uniformly distributed throughout the the Galactic disk. 
Early calculations by Hills (1972) indicated that the ion- 
izing radiation from these hot pre-white dwarf stars could 
significantly infiuence the interstellar medium. Since 
then, there has been progress in understanding the late 
stages of evolution of low mass stars, but little additional 
work has been carried out on their collective, large-scale 
infiuence on the interstellar medium. 



E. Is missing atomic data important? 



C. How much ionizing radiation escapes the galaxy? 

There have been several theoretical calculations of 
how much ionizing radiation escapes from galaxies (e.g., 
Ciardi et d] [2Q5^ [Clarke and Oeyl [2Q^| [Dove and 



Shullj |1994| IRicotti and Shull[ |2000| [Wood and Loeb| 
2000). The models show that if the interstellar medium 



is smooth the fraction of escaping photons is small. 
Clumpy density structures allow for larger escape frac- 
tions through low density paths in the interstellar 
medium (see Fig. [16] above). However, there are few ob- 
servations to test these models and determine how much 
ionizing radiation actually escapes. If interface radiation 
or shocks ( [Bland-Hawthorn et a/. , 2007) are not a ma- 
jor source of ionization, then one promising avenue is to 
study the ionized surfaces of distant high velocity clouds 
surrounding the Galaxy ( . Tufte et a/.||1998||2002D [Bland- 
Hawthorn and Maloney] (|1999D; [Putman al] d^QQSP; 



Tufte et al.\ (1998) have all estimated that a few per- 



cent of the ionizing luminosity from the Galaxy would 
be required to explain the Ha emission from high veloc- 
ity clouds located 
[2004) 



10 0^ pc above the disk (jWakker 
As discussed in 



jVj the spectrum of the escap- 
ing radiation (e.g., from future observations of He I/Ha 
from intermediate and high velocity clouds) could reveal 
whether photons escape the galaxy through very low den- 
sity channels or by filtering through density bounded ion- 
ized regions. 



Presently photoionization models are unable to make 
accurate predictions of the [S ll] emission from the DIG 
because the dielectronic recombination rates for sulfur 
are unknown (see discusion in Ali et aJ] [l991). The 
determination of these rates is important because they 
impact the inferred S+ and S++ abundances in the DIG, 
and, because sulfur is an important coolant, influencing 
the predicted strengths of other lines. While efforts are 
underway to calculate these dielectronic recombination 
rates, there may also be a handle on the rates obser- 
vationally via modeling the emission line strengths for 
large, low surface brightness H ll regions. 



F. What insights will new global models provide? 

A next step in modeling would be to test glo bal dynam- 
ical models of the interstellar med ium (e.g., de Avillez 
and Berry, 2001 Kowal et al [2007 ) to see whether their 
density structures can allow for O star radiation to pro- 
duce the observed ionization (Ha emission) and temper- 
ature structure (line ratios) of the DIG. This will be a 
formidable task, combining large-scale 3D dynamical and 
radiation transfer simulations. A unique solution for the 
structure of the interstellar medium may not even be 
possible without also incorporating observations of the 
other gas phases. However, combining dynamical and 
photoionization models would provide observational sig- 
natures that could be searched for. Perhaps progress 
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could be made in testing various scenarios — ruling out 
classes of models, determining what conditions (density 
and dynamics) are required to allow radiation to escape 
to the halo, and determining which models best fit the 
observed distribution and kinematics of the Ha over the 
sky. 

While much progress has been made in understand- 
ing the nature of the warm ionized medium and the 
basic physical processes occurring within it, the future 
promises to be a very exciting time, when the advances 
in computing ability are combined with high spatial and 
spectral resolution observations of this gas and other ma- 
jor phases of the interstellar medium. There is much still 
to be learned. 
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A list of acronyms and terms 

Chimney: a superbubble that extends from the galactic 
midplane, where it has been energized by super- 
novae, into the halo. Chimneys can be a conduit 
for hot gas and hydrogen ionizing radiation. 

DIG: diffuse ionized gas characterized by temperatures 
~ 10^ K and densities ^ 10~^ cm~^ that occupies 
the disk and halo of many spiral galaxies. In our 



Galaxy the DIG is often referred to as the warm 
ionized medium or WIM. 

DM: dispersion measure; the free electron density inte- 
grated along the line of sight to a pulsar; rie ds. 
Because the free electron density Ug ~ nH+ in the 
interstellar medium, DM is essentially the column 
density of the H+. 

EM: emission measure; the product of the proton den- 
sity times the electron density integrated along the 
line of sight. In regions of nearly fully ionized hy- 
drogen, emission measure is essentially J n^ ds. 

EPIC: European Photon Imaging Cameras on the 
XMM-Newton's X-ray telescopes. 

FIR: far infrared (wavelengths ~ 10^ /im). 

Forbidden line: an emission line produced by the decay 
of an excited metastable state in an atom or ion, 
denoted by square brackets around the ion's sym- 
bol. In the low density WIM/DIG, the collisional 
de-excitation time scales are much longer than the 
radiative decay time scales, and therefore the in- 
tensity of a forbidden line measures the excitation 
rate of the metastable state due to collisions by the 
thermal electrons in the plasma. 

the Galaxy: with a capital 'G', our galaxy, the Milky 
Way. 

Galactic latitude: The angular distance above the 
Galactic equator, the midplane of the Milky Way. 

Ha: the emission at 6563 A produced by the hydrogen 
Balmer-alpha (n = 3 ^ 2) transition following the 
recombination of ionized hydrogen. 

H I region: a cloud composed primarily of neutral hy- 
drogen atoms. 

H II region: a discrete region of photoionized hydro- 
gen associated with a hot star; a classical "emis- 
sion nebula" or "Stromgren Sphere" , as opposed to 
the more wide-spread, lower density H+ (the WIM, 
DIG) not clearly associated with a single, discrete 
source of ionization. 

H+: ionized hydrogen; in this paper ionized hydro- 
gen associated with the wide-spread, low density 
WIM/DIG. 

HST: the Hubble Space Telescope 

HVC: high velocity cloud; a neutral hydrogen cloud, 
usually located far from the Galactic disk and not 
partaking in the rotation of the disk. 

Ionization parameter: the ratio of the ionizing photon 
density to the electron density, which along with 
the spectrum, determines the population of ioniza- 
tion states in photoionized gas. 



25 



LIC: the very local interstellar cloud within a few parsecs 
of the sun. 

LyC: Lyman continuum; energies above the Lyman 
limit, the ionization potential of hydrogen (13.6 
eV). 

Magellanic Stream: an extended complex of neutral 
atomic hydrogen gas associated with a pair of 
satellite galaxies, the Large and Small Magellanic 
Clouds, orbiting the Galaxy. 

O and B stars: hot stars that emit Lyman continuum 
(LyC) photons. Massive O stars are located in ac- 
tive star forming regions, emit a large fraction of 
their luminosity in the LyC, and are believed to be 
the primary source of ionization for the WIM/DIG. 

PAH: polycyclic aromatic hydrocarbon; a large inter- 
stellar molecule that can inject heat into the inter- 
stellar via its ionization by ultraviolet photons from 
stars. 

R: rayleigh; a unit of surface brightness for emission 
lines, equal to 10^/47r photons s~^ cm~^ sr~^. An 
ionized region with an emission measure EM = 2.25 
cm~^ pc and a temperature of 8000 K has an Ha 
surface brightness of 1 R. 

SFR: star formation rate. 

Spiral galaxy: a galaxy that has a flat, rotating disk of 
stars, gas, and dust. 

Superbubble: a large cavity of hot (^ 10^ K), ionized gas 
created by the combined kinetic energy of multiple 
supernovae occurring within an active star forma- 
tion region. 

Tg; Tf. electron temperature and ion temperature, re- 
spectively. In general, Tg = Ti in the interstellar 
medium. 

TML: turbulent mixing layer; the transition region be- 
tween two adjacent parts of the interstellar medium 
that have very different temperatures and are mov- 
ing with respect to each other. 

VLA: The Very Large Array radio synthesis telescope 
near Socorro, New Mexico. 

WHAM: the Wisconsin H-Alpha Mapper; a remotely 
controlled observatory dedicated to the detection 
and study of faint emission lines from the interstel- 
lar medium of the Galaxy. 

WIM: the wide-spread, warm, ionized medium in the 
Galaxy characterized by temperatures ~ 10^ K and 
densities ^ 10~^ cm~^; sometimes also called the 
"Reynolds Layer". 

WNM: wide-spread neutral atomic hydrogen charac- 
terized by temperatures ~ 10^ K and densities 
- 10"^ cm-3. 



XMM- Newton: the Multi-Mirror Mission orbiting X-ray 
observatory of the European Space Agency. 

z: the perpendicular distance from the midplane of a 
spiral galaxy. 
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